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Abstract 

We provide a representation in terms of certain canonical functions for a sequence 
of polynomials orthogonal with respect to a weight that is strictly positive and ana- 
lytic on the unit circle. These formulas yield a complete asymptotic expansion for these 
polynomials, valid uniformly in the whole complex plane. As a consequence, we obtain 
some results about the distribution of zeros of these polynomials. The main technique 
is the steepest descent analysis of Deift and Zhou, based on the matrix Riemann-Hilbert 
characterization proposed by Fokas, Its and Kitaev. 

1 Introduction and background 

For r > 0, let B r be the open disc in C with center at z = and radius r, {z G C : \z\ < r}, 
and let T r be its boundary {z £ C : \z\ = r}. An integrable non-negative function w defined 
on Ti is called a weight if 

w(z) \dz\ > . 



For each weight w there exists a unique sequence of polynomials (p n (called Szego polynomials), 
orthonormal with respect to w, satisfying (f n (z) = K n z n + lower degree terms, K n > 0, and 



(p n (z)ip m (z) w(z)\dz\ = 5 mn ■ (1) 
Ti 
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We denote by $> n (z) = ip n (z)/K n the corresponding monic orthogonal polynomials. It is well 
known that they satisfy the Szego recurrence 

$ n +i(z) = z§ n (z) -c^$* n {z), <f> Q (z) = l, 

where we use the standard notation &* n (z) = z n $> n (l/~z). The parameters a n = — <£ n+ i(0) 
are called Verblunsky coefficients (also reflection coefficients or Schur parameters) and satisfy 
a n S Ox for n = 0, 1, 2, . . . (see Sim04a^ for details). 

The asymptotics of the sequence {<fn( z )} when n — > oo is well known for \z\ > 1 (the 
so called outer or exterior asymptotics) under rather general assumptions on the weight. 
However, the behavior of these polynomials in Bi is, in general, much more complicated and 
less studied. First, it is known that all the zeros of each <3?„ lie in U)i (see e.g. |GS841 p. 50-55], 
|Sze751 Ch. XII, theorem 12.1.1], |Ger61j or [Sim04a ). Moreover, they can be everywhere 
dense in the unit disc (the so-called Turan's conjecture), as shown in |AV88j . A recent result 
of Simon and Totik |ST04| proves even the existence of a "universal" measure of orthogonality 
on Ti for which the zero distributions of the corresponding {§ n (z)} approximate (in a weak 
sense) all the possible measures on the closed disc Bi. 

However, under stronger assumptions on the weight w we can say much more. In what 
follows we will consider only the strictly positive weights w on the unit circle that can be 
extended to an open set containing Ti as a nowhere vanishing holomorphic function. We 
will say briefly that such w is a strictly positive analytic weight on Ti. All the forthcoming 
statements are valid under this assumption. In order to shed light on this definition we 
point out that w is a non- vanishing holomorphic function in an annulus r < \z\ < 1/r (with 
< r < 1) and strictly positive on the unit circle Ti, if and only if, for any nonzero a£l 
the following representation is valid on Ti: 

w(z) = \W(z)\ a , zeTx, (2) 

where W is holomorphic and nowhere vanishing in the annulus. In particular, for a = 2, the 
following representation is useful: 

w (z) = W{z)W(z) = W{z)W{l/z) , z € Ti . (3) 

Obviously, here W is not uniquely defined. Normalized representations of this form arise 
from the so-called Szego function of w. Indeed, since w is strictly positive on Ti, its index 
(winding number about the origin) is 0, and we may define the Szego function (see e.g. |Sze751 
Ch. X, §10.2]): 

D(z) d = f exp (i- fj log W (e*) £±± tf) . (4) 

This function is piecewise analytic and non- vanishing, defined for \z\ ^ 1, and we will denote 
by .Di and D e its values for \z\ < 1 and \z\ > 1, respectively. It is easy to verify that 



Furthermore, if w satisfies the assumptions above, both D\ and D c admit a holomorphic 
extension across Ti, and maintaining the same notation for these analytic continuations we 
have 

AO) ( , f . 

D^zj = W{Z) ■ (6) 

This formula gives an analytic extension of w from Ti to an open domain containing the 
unit circle and, up to a normalization constant, provides the "standard" representation ©. 
Observe that by ©, w{z) = \D e (z)\~ 2 for z G T x . 
Let us define 

p = inf{0 < r < 1 : D c (z) is holomorphic in \z\ > r} . (7) 

The class of weights that we consider in this article is characterized by the fact that p < 1. 
Due to the symmetry ©, 1/p is the radius of convergence of the Taylor series for 1/Di about 
z = (hence, the constant 1/p coincides with Ftp-i in the notation of ShnO 1;\ . §7.1]). Taking 
into account ©, we may equivalently define p as 

p = inf{0 < r < 1 : l/w(z) is holomorphic in r < \z\ < 1/r} . (8) 

Again by ©, both circles T p and T 1 / p contain singularities of 1/w. 

Important results concerning the class of strictly positive analytic weights on Ti were 
obtained in |NT89j . It was shown there (see also |Sim04al §7.1]) that 



P 



lim |a n | 1/n . (9) 



In particular, strictly positive analytic weights on Ti are characterized by an exponential de- 
cay of the corresponding Verblunsky coefficients. Furthermore, the Szego asymptotic formula 

lim < (z) = — , zeB 1 , (10) 

can be continued analytically through the unit circle Ti and is valid locally uniformly in Bi/p. 
It shows that the number of zeros of {</? n } outside of D p remains uniformly bounded and these 
zeros are attracted by the zeros of D e in p < \z\ < 1 {Nevai-Totik points in the terminology 
of B. Simon SimO^Bj). Mhaskar and Saff MS90 , using potential theory arguments, gave a 
relevant complement to these results about zeros showing that for any subsequence {n^} C N 
satisfying 

p = lim|a„J 1/nfc , 

the zeros of {ip Hk+ i} distribute asymptotically uniformly in the weak-star sense on T p (see 
also |Pak87p . 

Further results on the asymptotics or zero behavior inside T p have been obtained so far 
only for some particular subclasses of strictly positive analytic weights on Ti. Besides the 
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classical explicit formula for the Bernstein-Szego polynomials (when w(z) = \l/q(z)\ 2 , q a 
fixed polynomial with zeros in C \ Ti), probably the first important result in this direction 
was obtained by Szabados |Sza79j . who considered the polynomial case w(z) = \p(z)\ 2 (where 
again p is a fixed polynomial with zeros in C \ Ti). Szabados gave a formula for the orthog- 
onal polynomials with respect to | 2 which, although containing unknown parameters, is 
sufficient for the asymptotic analysis inside T p . 

By extending Darboux's formula to Szego orthogonal polynomials, Ismail and Ruedemann 
|IR92j (see also |GM91j ) obtain matrix expressions for the 9? n 's for the case of rational weights 
w(z) = \p(z) I 'q(z)\ 2 , where the size of the matrix depends on the degrees of p and q. They 
applied their result to the asymptotic analysis in the case when T p contains a single zero of 
w. The same rational case is also the object of study in a less known paper il'akcT. , where 
the author considers the weak-star zero distribution for the <^ n 's. 

Another class of weights for which p < 1 corresponds to the Rogers-Szego polynomials, 
studied for instance, in 1MGH90] . where even the interlacing properties of their zeros have 
been established (see also |LS87j and |Sim04a| Example 1.6.5]). The reader interested in 
many more explicit examples of orthogonal polynomials on Ti is referred to |Sim04al Ch. I 
and VIII]. 

The main goal of this paper is to provide a complete asymptotic expansion for the orthog- 
onal polynomials, valid uniformly in the whole complex plane, under the assumption that w 
is a strictly positive analytic weight on Ti. As a consequence, we obtain results about the 
distribution of zeros of these orthogonal polynomials, under a variety of different assumptions 
about the weight function. These results recover, in particular cases, some of the results men- 
tioned above. The main technique is the steepest descent analysis of Deift and Zhou (see e.g. 
Dei99aJ), based on the matrix Riemann-Hilbert characterization proposed by Fokas, Its and 
Kitaev }FIK92| . This approach more generally allows one to handle analytic complex valued 
weights w that are nonvanishing on Ti provided the corresponding orthonormal polynomials 
<p n (z) exist with K n 7^ for each n. Furthermore, the asymptotic analysis is equally possible 
(and technically not much more involved) for a sequence of varying weights (when w depends 
on the degree n of the polynomial). This was carried out previously in BDJ99 for a specific 
weight relevant to the asymptotics of the distribution of the length of the longest increasing 
sequences in random permutations. A basic discussion of orthogonal polynomials on the unit 
circle, in the context of an example of integrable operators, and their connection with the 
Riemann-Hilbert problems is contained in [Dei99b . 

Until recently the analyticity of w has been a necessary condition for the applicability 
of the Riemann-Hilbert asymptotic analysis. However, McLaughlin and Miller MM04 have 
extended the approach via the solution of a d-bar problem. Their method allows one to 
handle piecewise smooth and strictly positive weights w on Ti; for this generality they pay 
the price of less detailed asymptotics. Also their results are qualitatively different: in the 
settings of MM04 essentially all the zeros approach Ti with a polynomial rate depending 
on the smoothness of w, and those that might stay inside are artifacts of the jumps of the 
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weight. 

Our asymptotic formulas (Theorem^) are stated in terms of a sequence of iterated Cauchy 
transforms of the function 

^)=A(^cW, (11) 

where D x and D e are respectively the inner and the outer Szego functions defined by (@J. By 
$7§, J 7 is holomorphic in p < \z\ < 1. Furthermore, by ©, 

T Q) = , for \z\ ^ 1 , |^(z)| = 1 on T x , (12) 

and l/^ 7 is holomorphic in 1 < \z\ < 1/p. 

The function l/T appears as a scattering matrix (or function) in |GC79| Formula (V.5)]. 
Geronimo and Case also give an asymptotic expression for <p n on Ti, where the Cauchy 
transform of the scattering functions appears implicitly (cf. formulas (HI. 7) and (V.8) in 
GC79 ). One of our goals here is to make precise the relevant role of the sequences of 
iterated Cauchy transforms related to T in the asymptotics of the orthogonal polynomials. 

The structure of the paper is as follows. In the next section we state the main results 
giving a representation of the polynomials <5 n , the leading coefficients K n and the Verblun- 
sky coefficients a n in terms of the series of some canonical functions. These formulas have 
an asymptotic nature that we exploit in Section El where we present a number of results 
concerning the asymptotic behavior of the Verblunsky coefficients, leading coefficients, and 
the "outer" asymptotic behavior of the orthogonal polynomials. We also present in Section 
121 general results describing the asymptotic behavior of the zeros, for a rather wide class of 
analytic weights. The detailed assumption is that the exterior function D e can be extended 
to the exterior of a domain {\z\ > p'} with finitely many poles. The proofs of these results 
are deferred to Sectional Finally, in SectionEJwe discuss some particular cases, including the 
case of a weight for which D e possesses an isolated essential singularity within the unit disk. 
The examples in Section [S] are all illustrated also with numerical experiments. We should 
point out that the zeros of $ n 's have been computed as the eigenvalues of the truncated ma- 
trix corresponding to the CMV representation of the multiplication operator (see [CMV03 
or Sim04a, Ch. IV] for details), using a symbolic algebra software with extended precision. 
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2 Canonical representation of the orthogonal polynomials 

Occasionally it is more natural to write our formulas in terms of a modified Szego function. 
Since w is of index 0, we may define the calligraphic functions 

„ , > dcf (I f log w(t)dt\ 

T>i(z) = exp — d> — , if \z\ < 1 , (13) 



_ 2m Jj 1 t — z 

V c (z) = exp — j—^ — J ' if |*| > 1 . (14) 

Both T>i and T> e are holomorphic and non-vanishing in the respective regions where they are 
defined and satisfy 

T ^m = ?M, H >i, 



Z J T 

where r is the Szego constant (mean value) given by 

f2w 



del 

T= D: 



k = 6XP ("S L '""fK^)*) — p(~ / T> log»W j) >0. (16) 
Furthermore (compare ©), 

P;(z)P c (z) = «,(*) . (17) 

This equation plus the normalization 2? e (oo) = 1 determines uniquely the nowhere vanishing 
functions H13|) - ()14j) . Finally, we have the following straightforward identities for \z\ ^ 1: 

Di(z) =TVi(z), if |^| < 1 , 

D c {z)=TV c ( Z y\ if|z|>l. 

In what follows, all the circles T Q , a > 0, are oriented counterclockwise. For a fixed value 
of a (which will be clear from the context), we set 

/+(z) d = f lim f(t), /_(z) d = f lim f(t), 

u — *z u — *z. 

\u\<a \u\~>a 

for any analytic function / for which these limits exist. We will also talk about the "+" and 
the "— " side of T a referring to its inner and outer boundary points, respectively. 

Choose an arbitrary r, p < r < 1, that we fix for what follows. We introduce the following 
operators acting on the space of holomorphic functions in C \ (T r U Ty r ) with continuous 
boundary values: 



2m J Tr w{t)(t-z) 2mr z J Tr t-z 

M c n {f){z) = —l /_(t) dt = —l f-(t)— — —. -dt, (19) 

nKJ,y> 2m J Ti/ / y! t^w{t){t-z) 2mJ Ti/r Jy) F(t)tn(t-z) 1 ; 



1 t tf . t n V({t) 1 I ,T(t)t n 
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where /_ denotes the exterior boundary values of the function / on the respective circles. 
Let us denote by 



A d = f ||^|| Tr =max|^(z)|; (20) 

zST r 



taking into account (|T2*|). also 

A 



1 

max 



T 1/r 



2ST 1/r \F{z) 



Then straightforward bounds show that 



\Ml(f)(z)\<^r n M^- z<£T r , 
1 1 r z \\z — r 



Ar 2 „ II/- 



\MUf)(z)\ <—r n "' " y , z £ T 1/P 
r ||z| — 1 /r | ' 



(21) 



Remark 1 It is obvious that for \z\ < p in (|18|) and for |z| > 1/p in ()19|1 we may take Ti 
as the contour of integration in these formulas. However, our choice makes inequalities (1211) 
evident. 



Using these operators we define the following sequence of functions: 

fk 0) = i, fP = K0), fW^MUfP), (22) 

and 

fV*» * KifW) , A 2k+2) d = f MM 2k+1) ) , ken. (23) 

Observe that each of these formulas defines in fact two analytic functions, one inside and one 
outside of the circle where the corresponding integral is taken. In order to avoid cumbersome 
notation we will not use the subindices i and e in this case, but the reader should bear in 
mind this fact. 



By CD, 



Ar 



\A+r |/(D| = |^(1)| <^ 



T \\Z\ 



A_2 „n A2„2n i 

|^l/r \A 2) \ = \MUA 1) )\<— TrT^T73l|/«||T 1/r < ' 



\z\*r => \fk 3) \ = \MM 2) )\<^jr^A\f^h r < 



X / r | IUn Pi/ r - I_ r || z | _l/ r | ' 
r A"V^ n i 



r 2 || z | _ r | M^n IHr - r 2 (I _ r ) 2 ||z| -r 
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and in general, for k €. N, 



, _. , if A; is odd, and \z\ 7^ r , 

t (1 — r ) \\z\ — r\ 

V, (24) 

A k r kn 1 



, if A; is even, and \z\ 7^ 1/r . 



Remark 2 There are a number of equivalent definitions for the functions ffl . For instance, 

-/y) ( z = <* 1 w dt = d) y -^dt = jy <* ^dt. 25 

V 7 2vri J Tr «;(<)(* -z) 2vri r 2 / Tr t - z 2vrir 2 J Tl t - z V 7 

We can rewrite the integrand directly in terms of the weight w recalling that on T l5 

T 2 V 71 "* Ai S-t J ' 

where the integral is understood as its principal value. If 

+00 

^(z) = ^ c fe z fc (26) 
fc=— 00 

is the Laurent expansion of T in the annulus p < \z\ < 1/p and V+ is the Riesz projection 
onto # 2 (Di), then for z £ D r , 

1 z n 
/W(z) = - — P + (z n .F(z)) = x (Laurent series of T truncated at z~ n ) , (27) 

(k) 

which establishes a connection of the fn s with Hankel operators. It also shows that for 
p <\z\ < r, 

n^oo z n T 

and convergence is uniform on each compact subset of this annulus. 
Finally, taking into account that 

it is easy to see that sequence {fn^} satisfies the following recurrence: 

/2iW=^ 1) -^. *eBi. (28) 
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Let 



^) = E/» M W. N^IA, and Z n (z) d = f £ /( 2fc+1 ) (ar) , N^r. (29) 

fc=0 fc=0 

Bounds 1)24(1 show that the series in the right hand sides of (|29(l converge absolutely and 
locally uniformly in the respective regions of definition if 

i iog(^F) 

Ar n < - - r n > i ^ . (30) 

r logr 

In this case there exists a constant C > depending on r and A only, such that 

C C 

\£n(z)\ < v;— i 7T7 ' for M 7^ V r > and \1n(z)\ < m r > for M / r • ( 31 ) 

||z| — l/r| ||z| — r| 

In the sequel we use C (eventually, with a subindex) to denote some irrelevant constants, 
different in each appearance, whose dependence or independence on the parameters will be 
stated explicitly. 

Now we may state one of the main results of the paper: 

Theorem 1 Let w be a strictly positive analytic weight on the unit circle T\, the constant p 
as defined in (Q) -(JEJ , and constant r with p < r < 1 fixed. Then with the notations introduced 
above, for every n S N satisfying (fHUl) the following formulas hold: 



r- 1 z n D c (z)£ n (z), if\z\>l/r; 



r-^z n D c (z)£ n (z)-^^-, ifr<\z\<l/r; (32) 
r1 n {z) 

— rTT~r ) if\z\<r. 



We emphasize that in each of the regions above £ n and I n have the meaning given in (|29l) . 
and are not obtained in general by analytic continuation from one domain to another. 

Remark 3 The canonical formulae expressed in (|32D appear to be valid for z bounded away 
from the circles of radius r and 1/r. However, Theorem ^ holds true for any r with p < r < 1, 
and so the formulae, in fact, are valid up to the boundaries, by analytic continuation. Indeed, 
by considering r with r < r < 1, the last formula in (|32j) is seen to remain valid by analytic 
continuation, for \z\ < f. Similarly, the first formula is valid by analytic continuation for 
\z\ > 1/r. For the middle formula in (|32j) . one considers r with p < f < r, which effectively 
extends the region where the middle formula is valid. Therefore, the formulae contained in 
121) are each in fact valid on the closure of the domains stated. 
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Remark 4 The second formula in (|32|) expresses the polynomial <£ n on Tj as a sum of 
"inner" and "outer" terms. An alternative two-term expression for <]?„ on Ti can be found 
in |GC791 Formula (III. 7)]. In the latter formula we can identify the inner and outer Szego 
functions, but the other factors are given implicitly in terms of the solutions of equations 
from the inverse scattering theory. 

Taking into account (|16ft . the following representation of the Verblunsky coefficients is 
straightforward: 

Corollary 1 With w as in Theorem^ we have for each n £ N, 

oo 
k=0 

Observe that under assumptions (|3(Jj) the series on the right hand side converges absolutely. 

The Riemann-Hilbert asymptotic analysis that we carry out in Section 0] allows us to find 
also a representation for the leading coefficients of the orthonormal polynomials. For this 
purpose we define the "shifted" sequence with respect to (|22 j) -(|23 |) : 

5 l 0)d = f i, ^ 1)d = f -M e „(i), gg^KW), 



and 



Inequalities ()21j) imply again that for k EN, 



j^2k r 2kn-\ 
.2k-l 



Theorem 2 With w as in Theorem^ and for n satisfying (|30j) . the leading coefficient k, 
of (p n has the following representation: 

2 00 
k=0 
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3 Asymptotic behavior 

By (|24l) it is immediate that for n satisfying ()3U|) and for N = 0,1,2, ... , 



N 



k=0 



< 



Cn 



\z\ — l/r\ 



„(2N+2)n 



\z\ / 1/r, 



(34) 



where the constant Cn depends only on r, N and A, but neither on n nor on z. Analogously, 

N 



k=0 



< 



Cn 



,(2AT+3)n 



(35) 



where Cn has a similar meaning as above. Here a similar observation concerning the validity 
of these bounds up to the boundary, made in Remark |31 applies. These bounds show that 
(|32fl allows us to obtain an asymptotic expression for {& n } of an arbitrarily high order. Let 
us discuss the consequences of truncating Z n and £ n in (|3*2"|) at their first terms. 

For instance, from Corollary ^ and Theorem [21 we get the following estimates of the 
Verblunsky coefficients a n and the leading coefficients K n : 

Corollary 2 Let w be a strictly positive analytic weight on the unit circle T±. With the 
notation introduced above and for each n G N, 



— / t n F{t)dt + (D(r 3n ) 
2m J Tr 



-c_ n _! + 0(r 3n ) , 



where c_ n _i is the corresponding Laurent coefficient of J 7 in (J^HJ)- Consequently, 

„3n\ 



\ < Ar n+1 + 0(r 



(36) 



Furthermore, 



2tt 



< 



A 2 r 



1 — r 2 



Another straightforward consequence of truncating I n and £ n at their first terms is the 
exterior asymptotics of {^n}: 

Proposition 1 Let w be a strictly positive analytic weight on the unit circle T\, and the 
constant p as defined in 0-©- Then 



lim 



n— >oo z 



(37) 
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locally uniformly for \z\ > p and convergence takes place with a geometric rate. In particular, 
for any compact set K in {z G C : \z\ > p and D e (z) ^ 0} there exists m = m(K) £ N such 
that for z S K , we have that 3> n (z) 7^ for all n > m. 

Furthermore, each zero Q of D e in p < \z\ < 1 attracts as many zeros of <F„ as its 
multiplicity, which tend to £ with a geometric rate. 

As mentioned in the Introduction, formula (JSTJ) (or the equivalent formula Q1UJI ) shows that 
in the analytic case the Szego's exterior asymptotics can be continued analytically through 
Ti up to the critical circle T p (cf. [NT89 or Sim04a, §7.1]). The rate of convergence is a 
consequence of (l3*2"|) . For instance, when r < \z\ < 1/r, by and (|3*H) - (|35|) . 



D e { Z ) 



DJz) 



(£n(z) ~ 1) 



a(*; 



■ln(z) 



< 



c 



(l/r-\z\)(\z\-r) 



where C depends on r and A only. 

The second statement of the Proposition is a direct application of this fact and of the 
Hurwitz's theorem. 

Regarding the zeros of $ n , we see that besides those matching the zeros of D e outside T p , 
the rest accumulates on D p . Nevai-Totik's theorem |NT89] states that in fact T p separates 
the regions of uniformly bounded and unbounded number of zeros of <3? n ; more precisely, 

p = inf {r > : 3 N(r) < oo so that Vn, § n {z) has < N{r) zeros in \z\ > r} 

(see also |Sim04al §7.1]). Furthermore, by l|3'2j). 



fW( z ) + h^(z) 



z G 



(38) 



and from ()34l) - (|35|) it follows that there exists a constant C, depending only on A and r, such 
that if p < r' < r, then 

C 



< 



,.3/f 



Z £ 



(39) 



As pointed out above, in a generic case the vast majority of zeros distributes uniformly at 
the critical circle T p , although some of them can remain inside, and by (|38|) - (|39D . stay close 

to the zeros of fn^ in D p . 

Let us make these statements more precise. For that purpose we must impose some 
additional conditions on our weight w (or function J-). 

Proposition 2 Assume that there exists < p' < p such that D e can be continued to the 
exterior of the circle T p /, as an analytic function whose only singularities are on the circle 
T p , and these are all isolated. Denote by a±, . . . ,a u the singularities (whose number is finite) 
of D c on T p , 

\ ai \ = ■■■ = \a u \ = p. (40) 
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Then for p < r' < r there exist constants < 5 = S(r') < 1 and C = C(r') < +oo such that 
for p' < \z\ < r' and n £ N, 



< C (p n 5 n + r in ) . 



(41) 



Furthermore, for every compact set K C 
C = C(K) < +oo such that for z G K, 



there exist constants < 5 = 5{K) < 1 and 



<$>n{z) 



1 y 



res .F(i) 

t=a fc \ t 



<C(p n 5 n + p 3n (l+£f n ) . 



(42) 



Remark 5 If T can be continued as an analytic function with a finite number of isolated 
singularities to whole disc H>i, then we may take 5 = in the right hand sides in (|41 |l — (|42 jl . 
Otherwise the right hand side in (|42jl may be replaced by an estimate of the form Cp n 5 n . 

Remark 6 This Proposition tells us that in general all the relevant information for the 
asymptotics of ^n's in D r comes from the singularities of the exterior Szego function D e on 
T p (that is, from the first singularities of D e we meet continuing it analytically inside the 
unit disc). 

Remark 7 In order to isolate the zeros of <3? n , one must be able to analyze the approximation 
to $ n afforded by (j41j) and (|42|) . For example, zero- free regions may be determined by (i) 
establishing zero- free regions for the approximation, and (ii) bounding $ n away from zero 
using the error estimates. Similarly, isolating the zeros can be done by first isolating the zeros 
of the approximation, and then using a Rouche' type argument for <£ n . 

Remark 8 It is interesting to note that Theorem ^ provides global asymptotics for <l? n , and 
as a consequence, Proposition [21 above provides an asymptotic description which is valid even 
in vicinities of the poles a^. 

Remark 9 It is an interesting challenge to find general assumptions concerning the measure 
of orthogonality that lead to an effective characterization of the asymptotic behavior of the 
zeros of § n . Simon Sim04b ( Sim04c] has placed assumptions on the Verblunsky coefficients 
which imply that all a^'s are simple poles. In this case the approximation afforded by (|41|) 
and (|42[) is completely explicit, and it is then straightforward (albeit technically complicated) 
to isolate and provide a complete asymptotic description for each and every zero of <& n . While 
we could carry out this extension of |Sim04b( ISim04c| , we rather choose to investigate a more 
general class of measures of orthogonality than considered in jSim04c| . and demonstrate how 
it is possible to obtain the same detailed level of asymptotic information concerning the zeros 
of <I> n . The first result in this direction is the following Theorem. 
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Definition 1 Let a, \a\ < 1, be a pole of a function / analytic in \z\ < 1. We denote by 
mult 2=a f(z) its multiplicity and say that a is a dominant pole of / if for any other singularity 
b of /, either \a\ > \b\ or \a\ = \b\, but then b is also a pole and 

mult f(z) > mult f(z) . 

z=a z =b 

In the sequel we use the following notation: for a £ C and e > 0, 

B e (a) = {z £ C : \z - a\ < e} . 

Theorem 3 Assume that there exists < p' < p such that D e can be continued to the 
exterior of the circle T p / ; as a meromorphic function whose only singularities are on the 
circle T p . Denote by ai,...,a u the poles (whose number is finite) of D e on T p , and assume 
that the dominant poles of D c are ai, . . . ,ag, t < u, and their multiplicity is m. 
Let e > 0. Then for p' < \z\ <r - e, z Ul =1 B £ (ak), and neN, 

*»(*) = ^ z n + ^7T E ( \) «r m+1 Di{ak)3ciak) + h n (z) , (43) 
r TDi(z) \m - 1J K a k - z 

where D e (ak) = lim 2 ^ afc D c (z)(z — ak) m , k = 1, ... ,£. There exist a constant < C < +oo 
independent of n and e, and a constant < 5 = 5(e) < 1, such that 



' C (p n 5 n + r 3n ) , ifm=l, 

K(z)\ < { C m _ 2 n (44) 
r n m z p n , if m> 2. 

g-m— 1 r ' J — 

Furthermore, for every compact set K C Bp there exists a constant C = C(K) < oo such 
that for z £ K , and ii£N, 



g 2ni(n— m+l)0fe 



e 



C5 n , if m = 1 , 
<<C .. >9 ( 45 ) 
n 

where ^ 

01 = 1, and 6» fe = — (arga fc - argai) , k = 2,...,£. (46) 

In particular, on every compact set K C D p , /or aW sufficiently large n polynomials <I> n can 
/iat>e ai mosi £ — 1 zeros, counting their multiplicities. 

Remark 10 In the spirit of Proposition [3 the asymptotics of $ n 's is governed by the dom- 
inant poles of J- only. However, the error term h n in Q45|) is small only away from all poles 
a k . 

Observe that this result is applicable to weights of the form w(z) = \R(z)S(z)\ 2 , z 6 Ti, 
where R is a rational function with at least one zero on T p (or one pole on ^i/p), and 5 is 
any function holomorphic and ^ in any annulus, containing {p < \z\ < 1/p}. 
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Remark 11 As mentioned in Remark |H] above, we may also obtain an asymptotic description 
of <& n in a vicinity of each of the poles a^, and describe the behavior of the zeros near these 
singular points. We leave this exercise to the interested reader. 

On the other hand, to investigate new phenomena, we consider in Section |S] the situation 
that there is a single isolated singularity of D c on the circle of radius p which is an essential 
singularity. For this case, we do carry out the detailed analysis of $ n near the location of the 
essential singularity. 

Corollary 3 Under assumptions of Theorem^ the Verblunsky coefficients satisfy 

^fn + l\ w _ m+1 ~ f0( P n S n ), ifm = l, 

an = ~Z^\ i) a k A(afc) D e {a k ) + < (47) 
£-^\m-ly K yO [n m p n ) , if m > 2 . 

Remark 12 The case of Verblunsky coefficients with the asymptotics (|47|) and m = 1 was 
analyzed in [Sim04c| . The subcase m = 1, u = 1 has been thoroughly studied in ^Si m04bj . 

Let us discuss now how Theorem |3] reveals the behavior of the zeros of the polynomials in 
D r . We distinguish two cases: the zeros approaching the critical circle T p and those staying 
inside D p . Under assumptions of Theorem 01 let us denote 

« y A(a fc )J>e(a fc ) e 2^ n - m +i)e k (48) 
ti ak ~ z 

9h defined in (|46|) . Numerical experiments show that that the zeros of approaching T p 

distribute radially uniformly, with gaps ("missing zeros") at the poles and at possible 
zeros of G n on T p . We will use the following notation: for a given analytic function / and 
for e > 0, 

Z(/) d = f {z: f(z) = 0}, 2 e (/) d = f |J B £ (z). (49) 

Also let 

B{e) = {z: p' + e<\z\<r-e, zi U u k=1 B £ (a k ) U Z £ (D C )} . 



Theorem 4 Let e > 0, and assume that the hypotheses of Theorem^ are satisfied. Put 

l/n 

,m — 1 



Vn= f p( n A ■ (50) 



Then there is a 5 > such that for all sufficiently large n £ N every pie slice of the form 

} 



z: v n (l-5/2) < \z\<v n (l + 5/2), - + ^ < arg (z) < - + ^ ) C g(g) \ Z e (G n ) 

n n 
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with appropriately chosen a E R, contains exactly k = k<i — k\ zeros of & n , z\ ,...,z k , 
satisfying 

and 

arg^-arg^^ + O^). (52) 




Figure 1: Zeros of $ n for n = 10, 25, 50, 75 (from left to right and from top to bottom) with 
w(z) = \(z — l/2) 10 (z — (j/2) 3 | 2 , z G Ti, where Q\ = exp(7ri\/2)- For comparison with the 
prediction (|5T|) we plot in each case with a dashed line the circle of radius ^ ( 1 + ^ log ( m ™ x ) ) • 



Remark 13 We may consider an example for the situation described in Theorem taking 
w(z) = \(z — l/2) 10 (z — Ci/2) 3 | 2 , z G Ti, with £i = exp(7ri\/2). The results of the numerical 
experiments appear in Fig. ^ Observe how the zeros exhibit the equidistribution pattern 
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close to the circle of radius 0.5 ^1 + ^ log { rr ^L^)J predicted by Theorem |31 and basically 
"ignoring" the non-dominant pole Ci/2 of T . An imaginative reader might observe though 
that the the influence of the non-dominant pole is not totally negligible: Fig. ^ shows a slight 
displacement of the zeros of & n towards the origin in a neighborhood of Ci/2. 

We can also be more precise about the accumulation set of zeros of ^n's on compact sets 
of D p . Let 

Z d = f D |J Z(<S> n ) . (53) 

k>ln>k 

Asymptotics 1)45(1 shows that the structure of Z n Bp depends on the relative positions of 
ai, ... ,ai on Tp. Without loss of generality we may assume that 9\ = 1, 9%, . . . , 9 V is the 
maximal subset of {6\, . . . , dg\ linearly independent over the rational numbers Q. Then there 
exist unique r k j G Q, k = 1, . . . , £, j = 1, . . . , v , such that 

V 

6k = i^ j r kj e j , k = l, ...,£, 

3=1 

(obviously, r^j = S k j, for 1 < k < v). Furthermore, we may always take r k i > 0, k = 1, . . . ,£. 
Denote 

rfci = — , 0<p k <q k eZ, k =!,...,£. 



Proposition 3 Let t G Z f]D p . Ifv = l (that is, if all 6 k G Q), then there exist < s k < qk< 
s^^TL, k = 1, . . . , £ such that 

k=l a k-t 

In this case Z n Bp is a finite set. 

If v > 2, i/ien additionally there exist X2, ■ ■ ■ , X v G R sttc/i i/iai 

A(afc)g e (afc) e 2 m ( Sfc / (?fe +E^ 2 r fcj x,) = Q f55) 

Remark 14 As it was observed by Szabados Sza79] (in the polynomial case), for v = 2 
equation (|55jl describes an algebraic curve of degree < £. In particular, if I = 2, Z n B p is 
a circular arc or a diameter in Bp (see Figure |2 left). For v > 2 equation (|55|) describes a 
two-dimensional domain bounded by algebraic curves (cf. Figure EJ right). Observe that we 
plot all the zeros of for n = 5, ... ,75, in order to reveal the structure of Z inside T^- 
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Figure 2: Zeros of $ n for n = 5, . . . , 75 with w(z) = \(z-l/2)(z - Ci/2)| 2 (left) and w{z) = 
|0 - 1/2)0 " Ci/2)0 - C2/2)| 2 (right), z G T 1; where (l = exp(vri^/2), C 2 = exp(m\/3). 



We can be more specific in the case of a rational weight, i.e. when w has the following 
form: 



w(z) = \r(z)\' 



1 



(56) 



where r is a rational function non vanishing on Ti; we may write 

r(z) = -(z)—(z), 



where p\ and q\ are polynomials without zeros in \z\ < 1, p e and q e are polynomials without 
zeros in \z\ > 1, and p c (0) ^ and g e (0) 7^ 0. Then, 



\r(z)\ 



Pi 1 \ Pe , \ 









|i?0) 



where i? has neither zeros nor poles in \z\ > 1; since in this identity we may multiply R by 
any power of z and any constant with absolute value equal to 1, without loss of generality 
we may assume that < R(oo) < 00. 



Definition 2 We say that w(z) = \R{z) 



1 is a normalized representation of the 



rational weight w if R is a rational function with neither zeros nor poles in \z\ > 1 and such 
that < R(oo) < 00. 
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If w(z) = \R(z)\ 2 , \z\ = 1, is a normalized representation of w, by rewriting this identity on 
Ti as 

w(z) = \R(z)\ 2 = R(z)Rjz) = R(z)R(l/z) , 
we immediately obtain that in this case, 

Di(z) = fl(l/2) , DtW-^, r = 1 ±- ) . (57) 

Furthermore, 

^)-¥r < 58 > 

i?(z) 

is a Blaschke product with its poles in \z\ < 1 coinciding with the zeros of R, and with its 
zeros in \z\ < 1 matching the poles of the same functions. 
Let ai,...,a£ be the zeros of R in \z\ < 1. Now by 

_ , , _ R(oo) R(oo) A (Ril/t) t n \ , , . , , 1 
R{z) R(l/z) ^ x t=a k \ R(t) t-zj 

where by (|H4jl-(|M5|l. |/i n (-2)| < C /{r — \z\) r 3ra , being C a constant depending on A and r only. 
In particular, if all the zeros of R in \z\ < 1 are simple, we get 

it(^) i?(l/z) R\a k ) a k - z 



Remark 15 As mentioned in the introduction, Szabados Sza79 obtained a formula similar 
to for the case of a polynomial weight; he also observed the equidistribution of the zeros 
approaching the critical circle T p . Rational weights with one dominant singularity have been 
discussed also by Ismail and Ruedemann in [JR92 . 

Remark 16 To contrast some of these results with a case of a more severe singularity, in 
Sectional we consider the weight of orthogonality 



w(t) = Wl(t) 



cxp 



a — t 



t 6 Ti 



(61) 



with < a < 1, as well as its inverse, 1/w. For w in the asymptotics of the Verblunsky 
coefficients is 



1 



a,. 



'2V5F 



3/4 



1 + 



n 



1 

T/2 



n 



00 , 
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(see Eq. p05[l ). where t + is the solution of the equation 



1 



1 



t 1 

+ 



t n + 1 \at — 1 t — a 



satisfying 



a + 



n + 1 \n 



n — > oo . 



where we take the positive square root. 

For the weight 1/w, with w defined in (|61|) . the asymptotics of the Verblunsky coefficients 

is 



(see Eq. (jl(J5|) ) , but t + is now the solution of the equation 



n — > oo , 



1 1 

+ 



t 1 

+ 



t n + 1 V at — 1 t — a 



satisfying 



f i = a + i 



n + 1 \ n 



where we take again the positive square root. 



n — > oo . 






0.5 



Figure 3: Zeros of $30 for weights w (left) and l/io (right), with w given in and a = 1/2. 
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Also the zero behavior in both cases is different. In order to entice the reader we show 
in Fig. |31 images of the zeros of (for n = 30 and a = 1/2), whose complete asymptotic 
description appears in Section |SJ 



4 Steepest descent analysis for orthogonal polynomials 



Observe that from the orthogonality conditions (0) it follows that 

if n {z)^w{z)\dz\ = - I ip n {z)z n - k - l ^-dz = Q for fc = 0, 1, . . . ,n - 1 . 



The reversed polynomials <p*(z) = z n ip ri {l/ z) also satisfy the following orthogonality condi- 
tions: 



<f* n _ 1 (z)z k ^-dz = i & ^ 1 {z)z k w{z)\dz\ 



0, fc = 0,l,...,n-2, 

i/Kn-i, k = n — 1 . 



Hence, 



Y(z) 



V 



1 £ $ n (t)w(t) dt \ 



2iri J Tl t n (t-z) 
K n-i I ¥*n-x{t)w{t) dt 



i J Tl t»(t-z) J 
is a unique solution of the following Riemann-Hilbert problem: Y is holomorphic in C \ Tj , 

'1 w(t)/t r 



Y + (t) = Y_(t) 



1 



z€T 1} 



and 



Inn He) °„ | I 



where I is the 2x2 identity matrix. Observe also that 

2 H(0) 
2tt ' 



(62) 



(see pBDJ99 or Dci99a] for more details on the Riemann-Hilbert characterization of orthog- 
onal polynomials). 
Let us define 

' °\ 7 11 

H(z) = M z n ) ' 1 |Z ' > ' (63) 



if \z\ < 1, 
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and put 

Then T is holomorphic in C \ Tj , 



T(z) = Y(z) H(z) . (64) 



T + (t)=T„(t) W t §) , teTx, 

and 

lim T(z) = / . 

2^00 

Next transformation is induced by the following factorization of the jump matrix for T: 

t n w(t)\ _ ( 1 OW w(t)\ ( 1 0\ 

t~ n ) ~ \l/{t n w(t)) l) \-l/w(t) ) \t n /w(t) l) ' 

Hence, with r, p < r < 1, fixed above we determine the following regions (Figure 0J: 

i~l = {z : \z\ < r}, fioo = {z : \z\ > 1/r}, 

+ = {z : r < \z\ < 1}, f2_ = {z : 1 < |z| < 1/r}. 
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We define now 
where 



U(z) d ^T(z)K(z), 



K{z) = { 



1 

z n /w{z) 1J 

1 o N 

l/{z n w{z)) 1, 



if z G $7o U Oqo, 

, if z e $7+, 

, ifzG^_. 



Then *7 is holomorphic in C \ (T r U Ti U T 1/r ) , 

lim U(z) = I , 



and 
where 



17+ (t) = U-(t) Ju(t), t G (T r U Ti U T 1/r ) , 



Ju(t) = < 



w(t)\ 
-l/w(t) ) 

1 o\ 

K t n /w(t) i y ' 

i o\ 

l/(t n w(t)) lj' 
Using the piece-wise analytic matrix 



if t G Ti, 
if t G T r , 
if t G 



1/r- 



2? e (^ 

■l/2?i(z) , 



£>e(«)/r \ 

t/A,(*) J ' 
( A(*)/r 



if 
if 



we make a new transformation, defining 

S(z) = U{z)N- 1 (z). 
Matrix S is holomorphic in C \ (T r U Ti U T 1/r ), lim^oo S(z) = I, and 

S+(t) = S-(t)J s (t), t G T r U Ti U T 1/r , 
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where 



J s (t) = N_J U N- 1 



1 -t n Vf{t)/w{t) 
1 

1 
[ \V*(t)/(t n w(t)) 



if t € Ti, 
, ifi€T n 

^ I , if t € T 1/r . 



This shows in particular that S is analytic across Ti. Using the definition of J 7 , the jump 
matrix function J$ can be rewritten as 



J s (t) =iV_J c/ iV+ 1 = < 



1 -t n ^(t)/r 2 
1 

1 0' 

I \vy(t»:F(t)) i 



if t € Ti, 
if t G T r , 



(69) 



if t € T 



l/r- 



Observe that 

J s (z) = J + 0(r n ) , Z £T r U T 1/r . 
Unravelling these transformations we have 

= S(z)N(z)K-\z)H(z) , 
where i?, K and iV are defined in (|63jl. (|66|) and l|67|l. respectively. 



(70) 
(71) 



Proof of Theorem In order to make formula (|71[) useful we need an asymptotic expression 
for S. It can be obtained if we observe that S solves the following Riemann-Hilbert problem: 
S + = S-Js on T r U T 1 / r , and S(z) = I + 0(1/ z) as z — > oo. The jump condition may be 
rewritten equivalently as 



(5 - /)+ - (S - I)- = S-(J 3 ~ I) on T r U T 1/r . 
Let us introduce the following notation: for / defined and continuous on T r U T-yi r 

f(t)dt 



(72) 



C(f)(z) = ^- i 



t - z 



is the Cauchy operator, and 



£-(/)(*) 



lim C(f){u) 



u — *z 

it6"-" side 



C+(f)(z) 



lim C(/)(«) 

iG"+" side 
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are the Cauchy projection operators on the "— " and "+"-sides of the corresponding curves, 
respectively. Furthermore, let 

V(f)(z) d = f C(f(J s - I))(z), V-(f)(z) d = f C_(/(J S - T))(*) ■ 

We can use the fact that Sokhotsky-Plemelj 's formula C + — C_ = I is valid also in 1? (cf. 
|Ste70j and |Dei99al §7.1]). 

Then, due to the normalization at infinity, (|72[) can be expressed, equivalently, as 

S - I = C(S-(J S - I)) = V(S-)=P(I)+V(S--I). (73) 

Considering the "— " boundary values in l|73j). we get that fx = S- — I satisfies fx = V-(I) + 
V-(li), or 

{I-V-)(ji)=V-{T). (74) 

This may be viewed as a singular integral equation for fx. Since C_ is a bounded linear 
operator in L 2 (T r U T]/ r ), there exists a constant < A < +oo, such that HP-/!! < 
for every / G L 2 (T r UT 1/r ), where we use the I? norm. In particular, using ([70 j) we get 

||P_(/)|| < A\\f(J s - < A|| J 5 - /||^ ll/H < Ar n \\f\\ , 

for n sufficiently large. For these n's, \\V- \\ < 1, and we can solve (|74|1 using Neumann series: 

oo oo 

M = {i-v-)-\r-V)) = E^(^-( J )) = E^-W' 

k=0 k=l 

or equivalently, 

oo 
fc=0 

Properly speaking, if we use the L 2 theory of Riemann-Hilbert problems, this series can be 
shown to converge in the I? sense. However, the series is, in fact, uniformly convergent. 
Indeed, straightforward contour deformations show that each successive integral in Vz.iT) 
taken over a circle of radius r may be deformed to a circle of radius r' with r' < r, and 
similarly each successive integral taken over a circle of radius \ jr may be deformed to a circle 
of radius 1/V with r' < r. Then arguments quite analogous to those leading to (|24jl may be 
used to estimate each term and the uniform convergence is evident. 

Replacing this series in (f73|) we get the full asymptotic expansion for 5: 

(oo \ oo oo 

fc=0 / fc=0 ' fc=0 
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where 

S<®=I, and S^=v[si k) 
Using the explicit expression for J$ we get 



2m 



(J s (t)-I)dt 



TrUT 1/r 



t - z 



k > 0. 



ti\z) 

9 { n\z) , 



where 



and operators M x n , M„ have been defined in (|18 [) -(|19 |) . Furthermore, 

S w (t)(J s (t)-I)dt 



S^(z)=V(S^)(z) = ^- 

zm ./T r UT 1/r 



t - Z 



1 



2™ fr 1/r \(9^)-(t) 







0\ dt 



V 2 e (t)/(t n w(t)) Ojt-z 



(fn ] )-(t)\ f0 -t n Vf(t)/w(t)\ dt 



2™ hr Wn ] )-(t) 







t - z 



T 

2m' jj 

r(2) 



l/r 



1 o\ (/P)-(t) * 1 



OJ F(t)t n t-z 2m r 2 



J) 





g { n\z) 



where 



d = f (Ml o A<«) (1) = MM ] ) , fi 2) = (M e n o Ml) (1) = M*(/«) 
Iterating this process, we get that 



/ 00 00 \ 



s* 



fc=l fc=0 

00 00 



\ k=0 



(2fc+l) 



(2fc) 
rt 

fc=l / 



where the matrix entries are defined recursively by 
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and 

It is usually the case in the asymptotic analysis of Riemann-Hilbert problems that the 
expansion for the error matrix is, at the end of the day, an asymptotic expansion, rather than 
being uniformly convergent. However, in the present setting, estimates (|24l) and (|33|) show 
that the series appearing in (|75|) are actually uniformly convergent (at least, for n satisfying 

Now we may replace ()75(l in (|71|) in order to find an expression for Y. We must do it 
independently in each region. 
In the domain floo, 

Hence, 

'z n D c (z)/r 

r/(z n D c (z)) 



Y(z) = S(z) 

and 

z n D c (z) 



^ n { z ) = Y ll {z) = S ll (z)z n D c (z)/r = ^^^ + J^jt k \z)j , \z\ > 1/r . (76) 
In the domain 0_, 

Hence, 

'z n DJz)/T 

w - i 

and for 1 < |z| < 1/r, 



y(z) " S(z) \-T/Di{z) r/{z n D c {z)) 



z n D (z) T 

* n {z) = Y n {z) = S 11 (z) ^ - S 12 (z) 



oo 



In the domain 



^)=( .Aw A( o )/T )' ^w^Uiw !)■ 



(77) 
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Hence, 

/ z n D e (z)/r AO)/r 



and we recover formula (|77l) . Hence, this asymptotic representation is still valid for r < \z\ < 
1. 

Finally, in the domain 0,q, 

"W-Uw A< o )/T )- KW = / ' " <2) = / ' 

Hence, 

^ = S WU W D,( o ,/T )' (78) 

and 

* n {z) = Y u (z) = -S 12 {z) ^j— = [f^f^Hz^j , \z\ < r. (79) 

This concludes the proof of Theorem ^ □ 

Proof of Theorem® The asymptotic expansion for the leading coefficient K n is a straight- 
forward consequence of formula (|62|) and the Riemann-Hilbert steepest descent analysis per- 
formed above. Indeed, by (|75)l and (|78j). 



T 522(2;) 

7)7~ I ' r 

k=l / 

and the statement follows from and (|62l) . □ 



Proof of Proposition® Assume first that \z\ < (p' + p)/2. Then t n J-{t)/{t — z) is a mero- 
morphic function of i in the annulus (p' + p)/2 < \t\ < r whose only poles are a±, . . . , a u . 
Hence, taking s = (3p + p')/4, we may deform the path of integration used in the definition 
of fn^ into the one depicted in Figure |3 left, obtaining that 



1 / Tit )t n , 1 -A / . , t n \ 1 / jf(* )t n , 

t-o <p r dt = — 9 7 res ^ (*) : - 7; — ; — 9 * , eft • 

it 2 J Tr t-z t * =afe V t-zj 2^r 2 J Ts i-z 



Since |z| < (// + p)/2, we have 

t n ^(t) 



t - z 



dt 



<47rmaxJ r (t)f n . (80) 

|t[=? 
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Hence, taking into account (|38 j) -(|39 |) we obtain (ji2*l) for \z\ < (p' + p)/2, with 5 = s/p. 
Obviously, as long as \z\ < s, (jlTj) is a consequence of (fl2"|) . 

Assume now that (p' + p)/2 < \z\ < r' . Then t n J-{t)/(t — z) is a meromorphic function 
of t in the annulus p' < \t\ < r whose poles are a±, . . . , a u and z. If z ^ {a%, . . . , a u }, the pole 
of t n ^(t) f(t — z) at t = z is simple, and 

res (V(i)T^) =z n F(z). 
t=z y t — z J 

Deforming the path of integration T r for fjp into the one depicted in Figure [51 right, with 
s = (3p' + p)/4, we obtain that 

/<■>(,) = -i/n.) - ^ E it -^{^Br dt - 

However, since \z\ > (p + p')/2, the same estimate Q8U|) for the last term in this identity 
remains valid. This proves (|41|). which remains valid after possible cancellations even for 
z G {ai, . . . , a u }. Finally, for any compact set K C Bp, we have that max z ^K \ z /p\ < 8 < lj 
and (021) follows from (jUJ). □ 

Proof of Theorem^ Observe first that if a, < \a\ < 1, is a pole of an analytic function / 
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in \z\ < 1, and in a neighborhood of z = a the Laurent expansion of / is 

(2 — a) m (z — a) m 1 



then 



/(t)t» = - a + a) n = /(*) f2 " °)' 

3=0 



and for n > m we get 

m— 1 



res (/(t)t») = £ Q a" >/ ; . = /_ m a""™ + 1 ^ (1 + K 



3=0 

where h n = for m = 1, and for m > 2 



(81) 



n |/-m| 

and constant Ci depends on a and m only. 

Thus, with assumptions of Theorem |3J for k = 1, . . . ,£, 

res = /_ m , fc (z) a"- m+1 f " ^ (1 + /^(z)) , (82) 

t=a k y t — z J \m — 1/ 

with 

tW _ , im nm-^r _ amw (83) 

t — 2! CLk — Z 

Furthermore, for z ^ B e (ak) we have that 

i , . C2 1 C2 1 

\h n h(z)\ < : : r < 7" , 

1 n ' U 71 ~ n Iz-afcl" 1 " 1 ~ n e™" 1 ' 

and constant C2 depends on p, r, A and m, but neither on z nor n. 

The multiplicity ml of a non-dominant pole at-, k > £, is strictly less than m; by (|81|) we 
obtain again that as long as 2; ^ B e (ak), 

res (V(t) = O ( p n n m '- x \ , m' < m , (84) 

t=a fc y t — Z J \ I 

uniformly in z. 
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Assumptions of Theorem |3] yield those of Proposition [3 Hence, plugging formulas l|82|l- 
into (|41 j) -l|42 |) . we obtain (|l3^) - (|4*H) . Estimate is a straightforward consequence of 
()43|) since we may assume now r 3 < p < r. 

For the last statement it is sufficient to observe that given a compact subset K in \z\ < p, 

by m, 

1 ( n 



a n-m+l \ m -l 



is a normal family in \z\ < p, and every limit of a convergent subsequence is a rational 
function with at most £ — 1 zeros on K. □ 

Proof of Theorem We carry out the proof in the spirit of |Sim04bj , but with an additional 
difficulty generated by the "floating" zeros of G n . 

Let e > 0; without loss of generality we may assume that a\ = p > 0. Then for z G B(e) 
we may rewrite (|4"3|l using the notation (flS]) : 



11 

m — 1 



*n(z) 



n 

m — 1 



^1 + 



G n (» 



p™- 1 .F(z) 



+ i? n (z), (85) 



where 



H n {z) 



D e (z) 



P 



n 

m — 1 



/i n (z) . 



According to there exists a constant C = C(e) such that 

C 



\HJz)\ < 



n 



z e B(e) . 



(86) 



Furthermore, {G n } is an equicontinuous family of rational functions of a bounded degree. In 
consequence, there exist constants q = q(e) and Q = Q{e) such that for all n G N, 



z£B(e)\Z £ (G n ) 0<q< 
and a constant <5 = 5(e) > such that 

z,(£B(e)\Z £ (G n ), \z-C\ <S 



G n {z) 



o m - 1 r{z) 



G n {z) 



<Q< +oo , 



G n (C) 



< 



p^Fiz) p m - l F{Q 
For v n defined in (JSUJ), there exists no G N such that for n > uq, 

v n <r-e, -<|, and (1 - 5/2) n < | < ^ < (1 + <5/2) n 
n 4 2 2 



(87) 
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Fix n > hq and assume that £, |£| = v n , is such that Bs(() C /3(e) \ Z £ (G n ). Denote 
G " (C) -se M and AJz) d = f Gn }^_ -se ia . 



Choose a neighborhood of £ = £ of the following form: 
/„ = I n (ki, k 2 ) 

„ / . , . , . . a + 2A;i7r , , a + 2/c27r 

zeC: v n (l - 6/2) < \z\ < v n (l + 5 2), — < arg (z) < — 

n n 

where integers k\ < k 2 satisfy £ 6 /„ and k 2 — k\ < dn/8. Then I n (ki, k 2 ) C B(e) \ Z e (G n ). 
By (jSTJ), |A n (z)| < s/4 on the boundary <9I n . In consequence, 

z€dl n \A n {z) + H n {z)\<\A n (z)\ + \H n {z)\<~, (89) 

where we have taken into account (|8fij) and (|88[). 

With the notation introduced above we rewrite as 

p- n ( U _ x ) = b n (z) + A n (z) + £T„(z) , (90) 

where 

For \z\ = v n (l-5/2), 

\b n (z)\>s-(l-5/2) n >s/2, 
where we have taken into account a condition in (|88j) . Analogously, when \z\ = v n (l + 8/2), 

\b n (z)\ > (l + 5/2) n -s>s/2, 

again by (|88|) . Finally, when arg(z) = (a + 2kn)/n, k £ Z, 

arg (( m -i) 1 ©")= a + 2to * fc-MI-U-i) "(7)" + s>s - 

In other words, we have that |6n( z )| > s/2 on z G (9/ n . Taking into account (|%9*|) . by Rouche's 
theorem, 6 n and 

/ \ -l 
t _„ / n 



D e (z) \m — 1 



32 



have the same number of zeros in I n (ki, k 2 ). But every pie slice I n (k, k + 1) of angle 27r/ra 
contains only one zero of b n , so in general <I> n has exactly k 2 — k\ zeros in every pie slice 
In(ki,k 2 ) C B(e) \ Z £ {G n ) of angle 2vr(A; 2 - kj)/n. 

On the other hand, if zq is a zero of <I> n , by (|§U)) . & n (zo) = — A n (zo) — H n (zQ), or 

m " x ) 1 (^) n = -*e--A n (z )-^o) 
so that 

- y'fwy_, +0 m + w=( ,(i + I, og f " )+©(! 

m — \J\pJ \n J \ n \m — 1 / \n 

Analogously, if 20,21 are the zeros of <l> n in two consecutive pie slices, I n (k — l,k) and 
I n (k, k + 1), then 

/ x 2vr / 1 

arg(zi) - arg(z ) = h C ~o 

n \ n , 

□ 

Proof of Proposition By (|4"5|) , 



z c p u z(g„) 



fe>ln>fe 

with G n defined in (|48|). and it is sufficient to describe all the possible limit points of {G n }. 
Observe that with the notation introduced before, 

G n (z) = > exp (27rz(n — m + 

fe=l flfc ~ Z 

^ Di(a k )D e (a k ) .-A \ 

= > exp Z7r? > rfcj (n — m + 1)9 j 

k=l ak ~ Z \ 3=1 ) 

The case v = 1 is trivial; assume v > 2. By Kronecker's theorem (also known as 
Kronecker-Weyl theorem, see e.g. |Cass571 Ch. Ill]), since 9\ = 1,9%, . . . ,9 V are rationally 
independent then for any real numbers X 2 , ■ ■ ■ , X v there exists a sub-sequence {rij} C N such 
that 

lim ^(K-m+l)^-^) = 1 jfe = 2, . . . , u , 

and (rij — m + l)p k /Qki k = 1, . . . ,£, have limits modulo Z. This shows that the set of limit 
points of {G n } is parameterized by < s k < q k , s k € Z, k = 1, . . . ,£ and X 2 , ... ,X V G R, 
and is given in the left hand side of (|55[) . Now the statement follows for u > 2. □ 
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5 Some examples 

We illustrate our results with some examples. Estimates (|3*%|) - ()39|) justify concentrating on 
the asymptotics of fn^ and of its zeros in D r . 

Consider first the simplest case of a polynomial weight inducing a single pole of D e inside 
TV Let w(z) = \z — a\ 2 , \z\ = 1, where without loss of generality we may assume a > 0, a ^ 1 
(the case a = is trivial). Then the normalized representation of w, according to Definition 
II is 



w(z) 



z — a 



Assume < a < 1; in this case 

D\(z) = 1 — az , D e (z 



if a < 1, and w(z) = |1 — az\ , if a > 1 . 

z (1 — az) 



z — a 



■ , r = 1 , and ^(z) 



z — a 



In consequence, with a < r < 1, 

fz n+1 (l-az)-a n+1 (l-a 2 ) 



a — z 



,n+l 



if |z| < r , 
if Izl > r . 



z — a 



Observe that for \z\ < r this is a polynomials of degree re, such that /n (a) 7^ 0, and whose 
zeros distribute asymptotically uniformly on Tui, leaving a gap at z = a. 
Furthermore, 

ff\z) = a 2n (l-a 2 ){l-az), z£B 1/n 

and in consequence, 



a — z 



Then by (jS2J), 



1 



z n+1 (l - az) - a n+1 (l - a 2 ) 



1 — az \ a — z 

z n+1 {l-az) 2 -a n+1 {l-a 2 ) 2 



+a 2n (l - a 2 



+ 0(a 5n ) ) , 



a — z 



which is valid locally uniformly in Bi a i +e , with e > 0. Obviously, we can get further terms in 
this asymptotic expansion. 
Analogously, 

1 — az ' 
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and 

n+l _ n n+l 

g M(z)=-a n (l-a 2 )- , zGB r . 

z — a 

Hence, using Theorem |2j we have the following asymptotics for the leading coefficient K n of 
orthonormal polynomials ip n : 



— (1 a 2n (l - a 2 ) + r 1 ")) , // - x . 



71 2ir 

Much richer is the case of a polynomial weight inducing two dominant poles of D e inside 
Ti. Consider 

w (z) = \(z - a){z - b)\ 2 , 0<a=|6|<l. 

Then 

D i (z) = (l-az)(l-bz), D c (z) = - y -, r = l, 

(z — a){z — b) 

and 

2 (i - Z )(i-M 

(z — a)(z — b) 

Let us denote b = ae tnd , and fix a < r < 1. Then we have for z6D r , 



f$\z) = -a n \r{z) (-) +a 



z\ n I - a 2 



a J 1 — e 



i-e 



1 ~ a e _ e iir(n+2)9 1 ~ a e 



a — z b — z 



(91) 



Hence, for \z\ < \a\, the zeros of fn are asymptotically close to the solutions of 

1 — n 2 p~ in6 , , 1 — n 2 P ilT 8 

1 Q e _ e in(n+2)e 1 a e 



a — z b — z 

or in other words, close to 



——7 — —rr , where A = . „ = „ . _ , \A \ = 1 . 92 



Let us consider some particular cases. First, 6 = y/2 ^ Q illustrates Remark 1141 In this 
case Z n Dui is parameterized by 



e"'" — a 
z = a — . a , , a G R , 

gi7ra Q—mO ' 



which is a straight line passing through the origin (cf. the numerical example in Fig. |^1 left). 
The rest of the zeros tends to the circle exhibiting the equidistribution pattern predicted 
by Theorem 0] (see Fig. EJ) . 
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Figure 6: Zeros of <& n for n = 10, 25, 50, 75 (from left to right and from top to bottom) with 
w(z) = \(z — 1/2) (2 — Ci/2)| 2 , z G Ti, where £1 = exp(7U\/2). The dashed line is the circle of 
radius 0.5. 

Next, if 9 = 1/2 (b = ia), we have by Proposition|3]that the zero accumulation set ZC\D\ a \ 
is finite. Namely, points (|92|) in this case have the following form: 

G D |o | if arg(A) G (-tt/4, 3tt/4) mod 2vr , 

G D| a | if arg(A) G (7r/4, 57r/4) mod 27r , 

G D| a | if arg(A) G (3vr/4, 7vr/4) mod 2vr , 

D| | if arg(A) G (5tt/4, 9vr/4) mod 2vr . 

Hence, ZnB| a | has at most 2 points, and unless 9 = ±tt/4 mod w, the sequence {fn^} 





n = 4k, k G N, 


> zo 


= ai 


i-A 
-I- A 


n = 


4k + 1, fcGN, = 


> Zl 


= ai 


-I- A 
-i-A 


n = 


4fc + 2, fc G N, = 


> Z2 


= ai 


-i-A 
I- A 


n = 


4fc + 3, k G N, = 


> Z3 


= ai 


I- A 
i-A & 
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-1 -0.5 0.5 1 



-1 -0.5 0.5 1 



Figure 7: Zeros of $ n for w(z) = \(z — l/2)(z — i/2)\ 2 , z G Ti, and n = 32, . . . , 35 (from left 
to right and from top to bottom). 



will have zeros in Bui for two consecutive indices n, and has no zeros inside T| a | for the next 

two. For instance, with a = 1/2, arg(^4) « 0.1567T, and zeros of {fn^} (and also zeros of 
{<!?„}) have limits in B| a | only for n = Ak and n = Ak + 3, k £ N, which are 

ai = i-(l+i), S2 = _!(! + ;), (93) 
lo fa 

(compare with the numerical experiment in Fig. [7J). 
On the other hand, by (j3K|) - (f33|l and by (JHU), 



Let us compare it with the results of BRLLS01 for 8 = 1/2. Denoting 



bj = lim _ ^' ty 7 n , > 3 = 1) • • • > 4 > 



$n(0) 
n=/"mod4 $„_l(0) 
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we have by (|94|). 



-i-A , 1-4 , i-il , -1-4 

Observe that |6 1 ...6 4 | 1/4 = |a|. Using the notation of BRLLS01 , we compute 



aS°(*) 



z + 6j abj + \ 



1 a + 6 



'i+i 



z& 





>j+2 



1 z + b j+2 

(where we take 65 = b\ and &6 = &2)> an d take 



A=|J{z: A^(z)=0}. 

For a = 1/2 it is straightforward to compute that AnB| a | = {si, S2}, where si, S2 are defined 
in ()93|) . I n other words, this case gives an example when the whole set A n Oui gives the 
accumulation points of zeros of {<& n } (cf. Remark 1 in [BRLLSO]]). 

Finally, we dwell very briefly on the case when D e has an essential singularity on T p 
analyzing two examples. 

Let first the weight function be given by 



w(t) 



cxp 



1 



a — t 



(95) 



with 0GB1. Without loss of generality we may assume < a < 1, so that p = a. Then 



Di(t) = exp 



t 



at -I 



, D c (t) = exp 



1 



t — a 



t = D (oo) 



1, 



and 



Tit) = exp 



t 1 

+ 



at — 1 i — a 

Again wee have to analyze the asymptotic behavior of the sequence of Cauchy transforms 

t n T{t) 1 r P »*n(t) 

f — z 2iri 



1 

2vri 



t - z 



dt , z £ 



(96) 



where we have used the notation 



= logt + - log^(z) = log* + - ( — 1— + — 

n n \at — It — a 
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Observe that we cannot apply the classical steepest descent method to the last integral in 
(|96j). since f„ depends on n. 

Let 7„ be a positively oriented closed Jordan curve encircling a and contained within B r , 
and eventually depending on n. Define 

If e^'Mi) 

'»M = W 7 „— *■ < 97 > 

Then 

^(i), j _ \-z n J r (z) + J n (z) , if z is in D r and outside j n , 
1 I n (z) , if z is inside 7„. 

For all sufficiently large n, function *$> n has four simple saddle points, t± and l/£±, 

0<t_<a<i + <l< l/t+ < 1/a < , 

such that 

t± = a± +0(l/n), n^oo. 



Function fy n is meromorphic in the cut disc Di \ (—1,0], and its level curves r = {t G 
C\(— oo, 0] : lm^ n (t) = Im , I' n (i + ) = 0} are trajectories of the quadratic differential *B' n (t)dt 
which has a double pole at t = with real leading coefficient, a pole of order 4 at a, and 
double zeros at t±. It is easy to see also that (0, 1) \ {a} C T. The typical local structure of 
the level set r at t = t+ is in Fig. EJ there are four arcs of T emanating from t+ with equal 
angles; two of these are (0, a) and (a, 1), another two should be vertical arcs passing through 
t± and ending at (0, 1). Since an analogous conclusion is valid in a neighborhood of we 
see that necessarily T contains a closed Jordan loop passing through t±, that we choose as 7„ 
in (|97|): see Fig. El right. Analogous arguments allow to determine the orthogonal trajectories 
r i = {teC \ (-oo,0] : Re* n (t) = Re* n (t+)}, see Fig. left. Since < *n(*+)> 

we have that Re^ n (t) < Re x I' n (t+) in the whole connected component of the complement of 

containing t-. In particular, Re x I / n (t) attains its strict maximum on 7„ at i+, so this is 
really the steepest descent curve for I n . 

Let us take 5 > and denote 



B n = B s/ ^{t + ) = {t G C : \t-t+\< S/y/K} , (3 n = dB n = {t G C : \t - t+\ = <5/v^}- 

(99) 

We assume n large enough, so that B n Cl r \ (—cxd, a]. Observe that t + is the only zero of 
* n (i) — fy n (t + ) in B n , and it is double: in a neighborhood of t = t + , 

* n (t) - * n (t+) = -^-^ (t - M 2 (i + o(t - , t -> t + . 
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a \ 












-*„(*+)) =0 ^ 











0.4 



0.5 



0.6 



0.7 



0.8 



Figure 8: Level curves Re(^f n (z) — i£ n (t+)) = and Jm{^ n (z)— \P n (i_j.)) = in a neighborhood 
of t± (for a = 1/2 and n = 30). 



By (HHD, *£(t+) = 2n 1 / 2 q - 3 / 2 (1 + o(l)). Condition y^(t) - tf„(t+) > for t > t + fixes a 
single- valued branch of \/^ n (t) — ^f n (t+) in B n . With this convention, we define 

Ut) = -in 1 / i y /* n (t)-V n (t+), (100) 

(where we take the positive square root). It is a conformal mapping of B n such that that 
Cn (K H -B n ) C £R, and Cn (7n H S n ) = (—d n , d n ) C R; moreover, the arc of 7 n n -B n in the 
upper half plane C + is mapped onto the positive semiaxis. Also 

n l/2 

Cn(t) = -i-g^ (t - t+) (1 + 0(t - t+)) , t^t+, 



and for sufficiently large n, 



S 2 

Cn(t) 2 >™, t£7n\^. (101) 
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0.8 - 
0.6 - 
0.4 - 
0.2 - 
- 
-0.2 - 
-0.4 - 
-0.6 - 
-0.8 - 
-1 - 



-0.5 



0.5 



Figure 9: Global structure of the level curves Re(^ n (z) — ^ n (t + )) = (left) and lm(^ n (z) 
*n(t+)) = for a = 1/2 and n = 30. 



For the error function 



i r+oo „-y/rit 2 

2vr« J_ 00 t - u 



1 n*n(t+)-\/«« 2 £ / 

2 V 



(102) 



i n 1 / 4 u 



we define F n (z) = f P n (( n (z)), z £ \ 7 n . Then F n has the following jump across 7 n PI B n : 
F ny+ {z) - F n _(z) = -e n *»^J e -v^Cn(^) 2 = _ e n*„(*) ^ 

matching the jump of across the same arc. This fact along with Sokhotsky's formulas 
show that 



Uz) 



I* n {z) + F n (z) - -L f ^-dt, for zeB n \ ln , 



F n (t) 

2m J Pn t- z 



dt. 



for z e c \ (b„ u 7 ;; 



where 



1 /• e n *™W 
F n (z) d ^-— J t ^—d1. ^ = ln \D„. 



t - z 



(103) 



(104) 
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Taking into account (|lUlj) . 



e nt> n (t+) / ^2 



Asymptotics for F n can be obtained using |AS72( Formula 7.1.23], 

e n*„(t+) a 3/4 / / 1 \\ 



also for z G £?„ \ 7„, 

< Cn" 1 / 4 e n ' I '' l ( i +) 



1 t KM* 



2iri J fat- z " dist(z, j3 n ) 

These estimates show that, roughly speaking, in the expressions (|10«-{j) for I n function F n 
dominates in B n , while its Cauchy transform along (3 n does it outside, and in consequence, 



In{z) 



t+-z 20FW V W /2 J) 



+ e 1 (z), if zeC\(£ n U 7 *), 



where 



f n (z) +£2(«), ifzS5 n \7 n , 



l£liz)l -6M(^) eXP {-^ n )> |£2(Z)I " dist(z,C„ 



This asymptotic expression has several corollaries. First, formulas (|98|) show that fn^ should 
have no zeros inside j n . All zeros lie outside j n and approach the circle T a with the clock 
pattern following asymptotically the solutions of the equation 

z\ n _ 1 /a\ 3/4 l ;F(t + ) 



t + J \nJ z — t + J-{z) 

Away from t + they should be close to the level curve Re(^ n (z) — ty n (t+)) = ^ log (^2^~J/i 
(compare with Fig. El left), although at the boundary j3 n of the disk B n they approach the 
level curve Re(¥ n (*) - ¥ n (M) = \ log (^f^) (see Fig. [TO] right). 

Furthermore, the Verblunsky coefficients a n for $ n exhibit the following asymptotic be- 
havior: 

an = -J—tlT(t + )(-) 3/4 (l + o(^)), n-00. (105) 



2v^ + v +/ \nJ V W /2 , 

A new feature is the fractional power of n in ()lU5f> . not appearing when the singularities of 
T on T p are only poles (cf. Corollary |3J) . 
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Figure 10: Zeros of <5 n for w given in Q95J) with a = 1/2 and n = 30, along with the level curves 
Re(¥ n (*)-tf n (i + )) = 0, Re(tf n (z)-¥ n (t + )) = ± log and Im(* n (z)-* n (t+)) = 

(left). Picture on the right shows that zeros of $30 in a neighborhood of the singular point a 
follow more closely the level curve Re(^ n (z) — \P n (i + )) = ^ log (77^^171)) as predicted. 



Consider now the inverse of the weight function given in (|95[) by setting 

2 



w(t) 



cxp 



1 



t — a 



, t€Ti 



(106) 



with < a < 1. Surprisingly we observe now a rather different behavior of the zeros of 
(cf. Fig. I13|) which is explained by the following analysis. For w in (|1U6|) . 



T{t) = exp 



t 1 

+ 



1 — at a — t 



and 

has saddle points 
satisfying 



= log* + - log^(z) =logt + - (-L- + — 
n n VI — at a — t 



t± —* a, l/t± — ► 1/a . 



f + = i_, t± = a ± i— = + 0(1 In) , n — ► 00. 
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Figure 11: Level curves T (thick lines) and T 1 - in a neighborhood of i+ for a = 1/2 and 
n = 30. 

Observe that in this case we have no dominant saddle point in a neighborhood of a. 

An analysis as for weight 1)95|) allows us to conclude that the local structure of the tra- 
jectories r = {ieC\ (-oo,0] : Im(* n (i) - V n (t+)) = 0} and r i = {teC \ (-oo,0] : 
Re(^I/ n (t) — \]/ n (i+)) = 0} is now as depicted in Figure ITT1 which yields in turn the global 
structure is as in Figure 1121 Now we take as j n in Q97JI the union of arcs of T lying in 
the domain where Ke{^ n (z) — \& n (i + )) < (see Fig. El right), except that we slightly 
depart from T in an arbitrary small neighborhood to the right of a, where the inequality 
Re(^ n (z) - * n (i+)) < is still valid. 

Along with the notation introduced in (|99[). (jlOOf) and l|lU2j) . we use and B* = {z : z E 
B n }, (3* = 8B* n , and 

F n (z) = P n ((n(z)), Z G fl„ \ 7 n, i 7 ^) = Pn(CnM), * G \ 7n- 
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Figure 12: Global structure of the level curves T (left) and F 1 - for a = 1/2 and n = 30. 



Following the arguments similar to those yielding (|103|) we obtain 



In(z) 



F n (z) + I*(z) 

F* n {z) + r n {z) 



F n (t) 



2ni J a t — z 



dt 



Kit) 



F n (t) 



2iri J a t — z 



F n (t) 



2m J a t — z 



dt 



dt - 

i r F*(t) 



2m la* t — z 

r'n 

1 



F*(t) 



2m la* t — z 

rn 



dt , for z G B n \ 7„ 
, for z G \ 7„ 



27ri /«» t — z 



for z£C \ (B n U Bl U 7^) , 
(107) 

where /* is as in (|104j) with 7* = j n \ (B n U5*). 

With contour 7„ described above we can perform the asymptotic analysis of the right 
hand sides in ()107j) : again F n and F* will dominate in B n and -B*, respectively, but the 
leading term outside these discs will be given by the sum of their Cauchy transforms. In 
particular, 



In{z) 



1 /a\3/4 /tWt+) t!f(t 



20F vn/ 



t_i_ — z 



+ 



t_ - z 



2ec\(s n uB:u 7 ;), 



and I n (z) has (asymptotically) at most one zero inside j n , originated by the cancelation of 
these Cauchy transforms. By (|98|) . the rest of the zeros of ffl lies outside j n and distributes 
asymptotically as the solutions of the equation 



z n T{z) 



1 /a\3/4 (t\T(t+) t n _T(t 



20F \n) 



t+ — z 



+ 



t_ - z 
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-1 -0.5 0.5 -1 -0.5 0.5 

Figure 13: Left: zeros of <!>„ for w given in with a = 1/2 and n = 30, contour j n , and 

the level curve Ke(^ n (z) — ^ n (t + )) = ^ log (2^^371) ( wrncn has two components). Right: 
same, but with n = 75. 

Hence, away from 7 n they should be close again to the level curve Re(^ n (z) — ty n (t + )) = 

n lo § (^S^) ( com P are with Fig.Hni). 

Finally, the Verblunsky coefficients a n of $ n satisfy 

Remark 17 Although we have considered here the case that D c possesses a single isolated 
singularity, it is clear that the same analysis carries through under the more general assump- 
tion that D c can be extended to a domain of the form {\z\ > p'} with p' < p such that D e 
possesses exactly one isolated essential singularity of the form (|96|) in this domain. 

Remark 18 Extensions in various directions of the considerations in this section are immi- 
nently possible; what is really required is (i) a domain {\z\ > p'} on which it is known that 
D c possesses a finite number of singularities and (ii) a description of the behavior of D c in 
a vicinity of each singularity. With this information, a complete asymptotic description of 
the function ' is within the reach of careful, and often creative, asymptotic analysis of 
integrals. 



46 



Acknowledgement 



The research of A.M.F. was supported, in part, by a research grant from the Ministry of 
Science and Technology (MCYT) of Spain, project code BFM2001-3878-C02, by Junta de 
Andaluci'a, Grupo de Investigation FQM229, and by Research Network "Network on Con- 
structive Complex Approximation (NeCCA)", INTAS 03-51-6637. A.M.F. acknowledges also 
the support of the Ministry of Education, Culture and Sports of Spain through the travel 
grant PR2003-0104, and the hospitality of the Department of Mathematics of the Vanderbilt 
University, where this work was started. 

The research of K.T.-R.M. was supported, in part, by the National Science Foundation 
under grants DMS-0451495 and DMS-0200749. 

The research of E.B.S. was supported, in part, by the U.S. National Science Foundation 
under grant DMS-0296026. 

Both A.M.F. and K.T.-R.M. acknowledge also a partial support of NATO Collaborative 
Linkage Grant "Orthogonal Polynomials: Theory, Applications and Generalizations", ref. 
PST.CLG.979738. 

We are grateful to Profs. P. Deift, J. S. Geronimo and P. D. Miller for very interesting 
discussions, as well as to Prof. B. Simon for his tireless activity in this field that stimulated 
our research as well. 



References 

[AS72] M. Abramowitz and I. A. Stegun. Handbook of Mathematical Functions. Dover 
Publ., New York, 1972. 

[AV88] M. P. Alfaro and L. Vigil. Solution of a problem of P. Turan on zeros of orthogonal 
polynomials on the unit circle. J. Approx. Theory, 53:195-197, 1988. 

[BDJ99] J. Baik, P. Deift, and K. Johansson. On the distribution of the length of the 
longest increasing subsequence of random permutations. J. Amer. Math. Soc, 
12(4):1119-1178, 1999. 

[BRLLS01] D. Barrios Rolama, G. Lopez Lagomasino, and E. B. Saff. Asymptotics of or- 
thogonal polynomials inside the unit circle and Szego-Pade approximants. In 
Proceedings of the Fifth International Symposium on Orthogonal Polynomials, 
Special Functions and their Applications (Patras, 1999), volume 133, pages 171— 
181, 2001. 

[CMV03] M. J. Cantero, L. Moral, and L. Velazquez. Five-diagonal matrices and zeros 
of orthogonal polynomials on the unit circle. Linear Algebra Appl, 362:29-56, 
2003. 



47 



[Cass57] J. W. S. Cassels. An introduction to Diophantine approximation. Cambridge 
Tracts in Mathematics and Mathematical Physics, No. 45. Cambridge University 
Press, New York, 1957. 

[Dei99a] P. A. Deift. Orthogonal polynomials and random matrices: a Riemann-Hilbert 
approach. New York University Courant Institute of Mathematical Sciences, New 
York, 1999. 

[Dei99b] P. A. Deift. Integrable operators. In Differential operators and spectral theory, 
69-84, Amer. Math. Soc. Transl. Ser. 2, 189, Amer. Math. Soc, Providence, RI, 
1999. 

[FIK92] A. S. Fokas, A. R. Its, and A. V. Kitaev. The isomonodromy approach to matrix 
models in 2D quantum gravity. Comm. Math. Phys., 147:395-430, 1992. 

[GC79] J. S. Geronimo and K. M. Case. Scattering theory and polynomials orthogonal 
on the unit circle. J. Math. Phys., 20(2):299-310, 1979. 

[Ger61] Ya. L. Geronimus. Orthogonal Polynomials: Estimates, Asymptotic Formulas, 
and Series of Polynomials Orthogonal on the Unit Circle and on an Interval. 
Consultants Bureau, New York, 1961. 

[GM91] E. Godoy and F. Marcellan. An analogue of the Christoffel formula for polynomial 
modification of a measure on the unit circle. Boll. Un. Mat. Ital. A, 5:1-12, 1991. 

[GS84] U. Grenander and G. Szego. Toeplitz forms and their applications. Chelsea 
Publishing Co., New York, second edition, 1984. 

[IR92] M. E. H. Ismail and R. W. Ruedemann. Relation between polynomials orthogonal 
on the unit circle with respect to different weights. J. Approx. Theory, 71(1):39- 
60, 1992. 

[LS87] D. S. Lubinsky and E. B. Saff. Convergence of Pade Approximants of Partial 
Theta Functions and the Rogers-Szego Polynomials. Constr. Approx., 3:331-361, 
1987. 

[MGH90] D. S. Mazel, J. S. Geronimo, and M. H. Hayes III. On the geometric sequences 
of reflection coefficients. IEEE Transactions on Acoustics, Speech, and Signal 
Processing, 38(10):1810-1812, October 1990. 

[MM04] K. T.-R. McLaughlin and P. D. Miller. The d steepest descent method and the 
asymptotic behavior of polynomials orthogonal on the unit circle with fixed and 
exponentially varying nonanalytic weights. Preprint arXiv:math.CA/0406484, 
2004. 



48 



[MS90] H. N. Mhaskar and E. B. Saff. On the distribution of zeros of polynomials 
orthogonal on the unit circle. J. Approx. Theory, 63:30-38, 1990. 

[NT89] P. Nevai and V. Totik. Orthogonal polynomials and their zeros. Acta Sci. Math. 
(Szeged), 53(l-2):99-104, 1989. 

[Pak87] L. Pakula. Asymptotic zero distribution of orthogonal polynomials in sinusoidal 
frequency estimation. IEEE Transactions on Information Theory, 33(4):569-576, 
July 1987. 

[Sim04a] B. Simon. Orthogonal Polynomials on the Unit Circle I and II. AMS Colloquium 
Series. American Mathematical Society, Providence, RI, 2005. 

[Sim04b] B. Simon. Fine structure of the zeros of orthogonal polynomials, I. A tale of two 
pictures, preprint, 2004. 

[Sim04c] B. Simon. Fine structure of the zeros of orthogonal polynomials, II. OPUC with 
competing exponential decay, preprint, 2004. 

[ST04] B. Simon and V. Totik. Limits of zeros of orthogonal polynomials on the circle, 
preprint, 2004. 

[Ste70] E. M. Stein. Singular Integrals and Differentiability Properties of Functions, vol- 
ume 30 of Princenton Mathematical Series. Princenton University Press, Prin- 
centon, 1970. 

[Sza79] J. Szabados. On some problems connected with polynomials orthogonal on the 
complex unit circle. Acta Math. Acad. Scient. Hung., 33(1-2):197-210, 1979. 

[Sze75] G. Szego. Orthogonal Polynomials, volume 23 of Amer. Math. Soc. Colloq. Publ. 
Amer. Math. Soc, Providence, RI, fourth edition, 1975. 



49 



